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ABSTRACT 

Results  from  localizing  Swallow  floats  using  range  data  collected 
during  the  8-9  July  1989  experiment  are  presented  herein.  As  part  of  the 
Downslope  Conversion  experiment,  the  Swallow  float  deployment  was 
conducted  near  34°  50'  N,  122°  20'  W,  about  150  km  west,  northwest  of 
Pt.  Arguello,  California.  Three  Swallow  floats  were  deployed  to  the 
ocean  bottom  and  9  freely  drifting  floats  were  ballasted  to  depths  starting 
at  600  m  and  spaced  every  400  m,  to  3800  m. 

Two  localization  methods,  least  squares  filter  and  Kalman  filter, 
were  applied  to  the  experiment  data  so  that  their  results  can  be  compared. 
Due  to  the  high  process  noise  experienced  by  the  freely  drifting  floats  dur¬ 
ing  the  experiment  and  the  quasi  vertical  line  array  deployment  geometry, 
the  two  filters  performs  comparably.  The  rms  position  error  is  estimated 
to  be  less  than  3.3  m  for  the  Kalman  filter  and  less  than  4.6  m  for  the  least 
squares  filter.  Both  filters  appear  to  be  capable  of  estimating  float  posi¬ 
tions  to  within  the  desired  accuracy  of  one-tenth  of  a  wavelength  at  the 
highest  frequency  of  interest  20  Hz  (7.5  m)  in  order  to  effectively  beam- 
form  the  VLF  acoustic  data. 
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1.  Introduction 


This  report  presents  the  results  from  localizing  Swallow  floats  using  range  data 
collected  during  the  July  1989  experiment  which  took  place  between  8  and  9  July  near 
34°  50'  N,  122°  20'  W,  about  150  km  west,  northwest  of  Pt.  ArgueUo,  California.  In  this 
experiment,  12  Swallow  floats  were  deployed;  9  were  freely  drifting  in  the  water  column 
in  a  quasi  vertical  line  array  geometry  with  a  vertical  float  separation  of  about  400  m, 
starting  at  about  600  m  deptli  to  about  3800  m;  and  3  were  tethered  to  the  ocean  bottom 
so  that  stable  absolute  position  fixes  could  be  obtained  for  the  freely  drifting  midwater 
floats.  Detailed  information  about  the  experiment  is  contained  in  the  experiment  trip 
repon  [1]. 

The  theoretical  basis  and  simulation  results  for  the  two  localization  methods, 
least  squares  filter  and  Kalman  filter,  can  be  found  in  [2,3].  A  comprehensive  procedure 
for  estimating  float  positions  using  data  from  the  April  1987  experiment  was  documented 
in  [4].  The  emphasis  here  is  on  the  use  of  the  filters  and  interpretation  of  the  results  using 
data  from  the  July  1989  experiment. 

This  report  is  organized  as  follows.  Following  this  introduction  section.  Section 
2  reviews  the  two  filters  from  the  performance  criteria  point  of  view.  The  desire  here  is 
to  provide  the  necessary  background  for  discussion  in  the  subsequent  sections.  Section  3 
describes  the  inputs  to  the  filters  while  Sections  4  and  5  discuss  approaches  to  filter  tun¬ 
ing  and  the  outputs  from  the  filters  respectively.  Lastly,  summary  of  the  localization 
results  is  given  in  Section  6. 
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2.  Loc..lization  Methods 

The  function  of  the  localization  filter  is  to  algorithmically  compute  float  posi¬ 
tions  from  noisy  range  measurements  in  a  statistical  sense  depending  on  the  optimality 
criterion  chosen.  In  this  section,  two  such  filters  developed  at  the  Marine  Physical 
Laboratory  (MPL)  are  reviewed  from  the  performance  criteria  point  of  view. 

2.1.  Generalized  Least  Squares  Filter 

The  generalized  (or  weighted)  least  squares  filter  [5,6]  involves  using  the 
current  set  of  range  measurements,  Z„,  which  are  related  to  the  set  of  float  positions  by 
the  expression: 


Zn  =  h{Xn)  +  Vn  (measurement model)  (2. 1 ) 

where  h  is  the  nonlinear  function  which  gives  the  ideal  (noiseless)  connection  between 
the  range  measurements  and  the  float  positions,  and  Vn  is  the  set  of  range  measurement 
errors.  These  errors  are  also  known  as  the  residuals  with  zero  means  and  second-order 
statistics  described  by; 


E 


=  Rn 


(measurement  statistics) 


(2.2) 


i.e.  the  errors  are  mutually  uncorrelated. 

The  least  squares  method  is  concerned  with  determining  the  most  probab’e  set 
of  float  positions  which  is  defined  as  the  set  of  float  positions  that  minimizes  the 
weighted  sum  of  the  squares  of  the  residuals: 


rmni. 


- 

ijpuze  Zn-hiXn)  ’’’  Wn 


Z,  -/i(X,)l 


where  W,  is  the  weighting  matrix  chosen  as  the  inverse  of  the  range  covariance  matrix 
Rn ,  and  the  weighted  least  squares  filter’s  position  estimate  is: 
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Xn 


HT  Rn  H  '*  Hr  /?,  Z„ 


(2.3) 


where: 


H  = 


dhjX^)  _  dh 
dX  dX 


X. 


(2.4) 


Thus,  the  least  squares  filter  estimates  current  float  positions  using  only  current 
set  of  measurements  and  an  estimate  of  the  measurement  error  statistics. 


2.2.  Kalman  Filter 

The  Kalman  filter  [3,6, 7, 8,9]  attempts  to  better  estimate  float  positions  by  tak¬ 
ing  advantage  of  the  knowledge  of  float  dynamics  and  incorporating  a  fading  memory  of 
past  data  into  the  estimator  structure.  In  addition  to  the  measurement  model  (2.1)  and 
statistics  (2.2),  Kalman  filter  also  uses  the  system  model  (float  dynamics),  the  system 
statistics,  and  the  initial  conditions.  The  system  model  is  given  by: 


X,  =  <t>  Xn-i  +  r  Xn-i  (system  model)  (2.5) 

where  d)  is  the  state  transition  matrix  which  relates  X„_i  to  X„,  r  the  matrix  which  relates 
the  accelerations  to  the  float  positions,  and  Xn-i  the  set  of  accelerations  with  zero  means 
and  second-order  statistics  described  by: 


E  X„  XI 


=  Qn  dm  (system  statistics) 


The  initial  conditions  are: 


(2.6) 


Xo  =  E(Xo) 


(2.7) 


Po  =  E{{Xo-Xo){Xo-Xo)r 
L 


(2.8) 


where  Xo  is  the  estimate  of  the  initial  float  positions,  and  Po  the  estimate  of  the  error  vari¬ 
ance  in  the  initial  float  positions. 
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The  operation  of  the  Kalman  filter  can  be  viewed  as  a  predictor-corrector  pro¬ 
cess  [10].  First,  suppose  we  have  available  to  us  the  previous  float  position  estimates  Xn-\ 
and  associated  position  covariance  matrix  P^-x: 


P.-i  =  E 


(2.9) 


and  we  would  like  to  obtain  the  best  estimate  of  the  current  float  positions  based  on  the 
previous  estimate.  We  are  in  the  "prediction  phase"  of  the  process.  The  system  model 
(2.5)  is  used  to  predict  the  float  positions  and  the  associated  position  covariance  matrix: 


X.(-)  =  <l>X,-i  (2.10) 

P,(-)  =  4>Pn-i<t>r  +  rx,  P-  (2.11) 

The  "minus"  is  introduced  here  (following  the  notation  in  reference  [6])  as  a  reminder 
that  this  is  our  best  estimate  prior  to  incorporating  the  measurement.  We  now  seek  to  use 
the  measurement  z,  to  improve  the  estimate  x,(-).  We  choose  a  linear  blending  of  the 
noisy  measurement  and  the  X,(-)  in  accordance  with  the  equation: 


X„  =  X, (-)  +  X,  [z,  -  /I  (X,  (-))]  (2.12) 

The  optimal  choice  of  X, ,  the  Kalman  gain,  is  to  minimize  a  scalar  sum  of  the  diagonal 
elements  of  the  position  covariance  matrix  p, : 


where: 


minimize 

K. 


' 

-  trace 

p. 

j 

Pn=E 


(X, -X,)(X, -X,)r 


(2.13) 


X,  is  obtained  as: 
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or  in  a  simpler  form: 


Kn  =  Pn (-)  HI  [//,  Pn  (-)  Hi  +  /?„ 


-1 


(2.14) 


where: 


Kn=PnHlRn^ 


(2.15) 


u  —  dh  (Xn  ) 
" — 


(2.16) 


|x.=x.(-) 

We  then  use  the  current  set  of  range  measurements  Z,  to  determine  the  innovations.  The 
innovation  is  defined  as: 


which  is  the  difference  between  the  actual  range  measurement  and  the  predicted  range 
measurement  (2.10).  Now  we  enter  the  "correction  phase"  of  the  process.  That  is,  we 
correct  or  update  the  predicted  positions  based  on  the  new  information  in  the  range  meas¬ 
urements  -  the  innovation.  The  innovation  is  weighted  by  the  Kalman  gain  Kn  to  correct 
the  predicted  positions  X,(-)  as  described  in  (2.12).  The  associated  position  covariance 
matrix  (2.1 1)  is  corrected  as  well: 


/  -KnHn 


Pn{-) 


(2.17) 


The  predictor-corrector  process  then  repeats  until  all  measurements  are  consumed. 

A  useful  interpretation  [6]  of  the  Kalman  gain  (2.15)  is  that  the  gain  is  "propor¬ 
tional"  to  the  uncertainty  in  the  position  estimates,  and  "inversely  proportional"  to  the 
range  measurement  noise.  In  other  words,  if  the  range  measurement  errors  are  large  and 
the  predicted  position  errors  are  small,  the  innovation  is  due  chiefly  to  the  noise  and  only 
small  change  in  the  predicted  positions  should  be  made.  On  the  other  hand,  small  meas¬ 
urement  noise  and  large  uncertainty  in  the  position  estimates  suggest  that  the  innovation 
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contains  considerable  information  about  errors  in  the  position  estimates.  Therefore,  the 
difference  between  the  actual  and  the  predicted  range  measurement  will  be  used  as  a 
basis  for  strong  corrections  to  the  position  estimates. 

In  short,  the  Kalman  filter  estimates  float  positions  by  weighting  a  current  set 
of  measurements  against  previous  position  estimates  propagated  forward  in  time  using  an 
equation  of  motion.  Thus,  the  Kalman  filter  uses  current  measurements  as  well  as  past 
measurements  in  a  fading  memory  fashion,  estimates  of  the  measurement  error  statistics 
and  acceleration  error  (process  noise)  statistics,  and  the  float  dynamics.  It  has  been 
shown  [3]  that  the  Kalman  filter  outperforms  the  least  squares  filter  in  the  presence  of 
high  measurement  noise  and  when  the  process  noise  is  low  because  of  its  ability  to  track 
float  motion  and  effectively  smooth  noisy  measurements. 
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3.  Inputs  to  localization  filter 

This  section  presents  the  inputs  to  both  filters  using  data  from  the  July  1989 
experiment.  Based  on  discussion  in  the  previous  section,  input  data  falls  under  four 
categories:  measurement  data,  measurement  statistics,  initial  estimates  and  process  noise. 

3.1.  Measurement  Data 

Range  measurements  used  by  the  filters  are  derived  quantities  computed  from  8 
kHz  pulse  travel  time  measurements  between  floats  (or  from  float  to  surface)  using  an 
estimate  of  the  sound  speed.  Thus,  measurement  data  includes  travel  time  estimates  and 
sound  speed  estimates. 

3.1.1.  Travel  Time  Estimates 

Measurement  of  pulse  travel  times  has  received  a  great  deal  of  attention  during 
the  process  of  the  Swallow  float  system  development  in  the  past  few  years.  References 
[1,11,12]  have  addressed  this  subject  in  detail.  In  brief,  each  float  is  equipped  with  an 
acoustic  transducer  which  transmits  and  receives  8  kHz  pulses.  Pulses  are  transmitted  by 
the  floats  in  a  preprogrammed  sequence.  A  different  float  transmits  every  45  seconds;  12 
floats  were  deployed  in  this  experiment  and  each  float  transmits  every  9  minutes.  A  given 
float  receives  pulses  transmitted  by  other  floats  as  well  as  its  own  arrivals  (surface 
echoes,  bottom  bounces,  or  mixtures  of  both).  Since  only  direct  path  pulses  between 
floats  and  surface  echo  pulses  from  one  float  to  itself  are  of  interest,  an  edge  detector  pro¬ 
gram  [11,3,4]  was  used  to  detect  and  extract  such  pulses  which  correspond  to  the  first 
arrivals  in  a  narrow  range  time  window. 

Surface  echo  travel  time  measurements  are  obtained  by  subtracting  outgoing 
pulse  transmit  times  from  pulse  arrival  times,  whereas  interfloat  travel  time  measure¬ 
ments  are  obtained  by  subtracting  pulse  transmit  time  according  to  the  transmitting  float 
from  pulse  arrival  times  at  the  receiving  float.  Since  two  time  bases  are  involved. 
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interfloat  travel  time  measurements  contain  a  bias  due  to  variation  in  the  float  clock  rates. 
It  has  been  shown  [12]  that  the  bias  can  be  reduced  significantly  by  averaging  reciprocal 
path  interfloat  travel  time  measurements.  Because  interfloat  measurements  were  not 
made  simultaneously,  travel  time  measurements  must  be  interpolated  by  a  factor  of  12 
before  reciprocal  path  measurements  can  be  averaged. 

Figures  3.1  through  3.12  show  the  leading  edge  of  each  float’s  detection  of  its 
own  transmitted  and  surface  echo  8  kHz  pulses.  The  vertical  axes  in  the  figures  have 
been  scaled  from  travel  time  to  depth,  using  half  of  1500  m/s  for  sound  speed.  The  line  of 
dots  at  approximately  0  depth  corresponds  to  the  detection  of  the  outgoing  8  kHz  pings 
and  the  second  line  of  dots  corresponds  to  the  detection  of  the  incoming  pulses  i.e.  sur¬ 
face  echo  pulses.  Note  that  the  lines  of  dots  terminate  prematurely  for  floats  4  and  5  in 
Figures  3.5  and  3.6,  respectively.  Roat  4  stopped  recording  data  4  hours  too  early  and 
float  5’s  automatic  ballast  release  time  was  set  incorrectly  [1].  Figures  3.13  through  3.78 
(except  Figures  3.13,  3.14,  and  3.24)  show  the  direct  path  pulse  detections  between 
floats.  Figures  3.13,  3.14  and  3.24  show  the  surface  reflected  pulse  detections  between 
the  bottomed  floats,  9,  10  and  11.  The  direct  path  pulses  between  the  bottomed  floats 
were  not  detected  because  the  sound  speed  gradient  bends  nearly  aU  of  the  sound  energy 
upward  over  the  6.3  km  distance  which  separates  the  floats.  Direct  path  travel  times  can 
be  calculated  from  surface  reflection  travel  times  using  equation  (3.1)  [4]: 


2 

surface  reflection 


-  T» 


(3.1) 


where  x,"  is  the  interfloat  travel  time,  t,;  and  x„  the  individual  float  surface  echo  travel 
times,  fa.*  the  harmonic  mean  sound  speed  for  a  vertical  path  from  the  surface  to  the  bot¬ 
tom,  and  c,,  the  sound  speed  at  the  depth  of  floats  i  and  j. 


Figures  3.79  through  3.87  contain  the  surface  echo  travel  time  measurements 
for  the  freely  drifting  floats.  Figures  3.88,  3.89  and  3.90  show  the  surface  echo  travel 
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time  measurements  for  the  bottomed  floats,  9,  10  and  11.  As  pointed  out  in  [4J,  the  bot¬ 
tomed  floats’  surface  measurements  like  these  are  expected  to  be  approximately  constant 
because  the  floats  themselves  are  approximately  stationary;  they  are  tethered  to  the  ocean 
bottom  by  3.05-meter  lines.  Noise  in  the  surface  echo  pulse  arrival  time  is  thought  to  be 
caused  by  scattering  of  the  pulse  at  the  rough,  moving  sea  surface  and  destructive 
interference  among  multiple  arrivals  at  the  receiver.  Since  filter  performance  can  be 
improved  by  substituting  a  constant  for  the  actual  measurements,  the  most  probable  value 
i.e.  the  mode  of  the  travel  time  measurements  (5.402  seconds  for  float  9,  5.422  ..cconds 
for  float  10,  and  5.404  seconds  for  float  1 1)  will  be  used  as  the  travel  time  estimates. 

Figures  3.91a,  3.92a,  and  3.93a  show  the  averaged  surface  bounce  path  travel 
times  between  the  bottomed  floats.  Like  the  bottom  float  surface  echo  travel  times,  the 
interfloat  travel  time  should  be  approximately  constant  because  the  floats  are  approxi¬ 
mately  stationary.  Figures  3.91b,  3.92b,  and  3.93b  show  the  direct  path  travel  time  esti¬ 
mates  between  the  bottomed  floats,  calculated  using  equation  (3.1)  and  the  actual  surface 
echo  travel  times,  whereas  Figures  3.91c,  3.92c,  and  3.93c  show  the  direct  path  travel 
time  estimates,  calculated  using  equation  (3.1)  and  the  constant  value  (mode'  estimates 
of  the  surface  echo  travel  times.  Because  the  underlying  direct  path  travel  time  must  be 
approximately  constant,  the  modes  of  Figures  3.91c,  3.92c,  and  3.93c  will  be  used  as  the 
direct  path  travel  time  estimates  (4.025  seconds  for  path  between  floats  9  and  10,  4.157 
seconds  for  path  between  floats  9  and  11,  and  4.132  seconds  for  path  between  floats  10 
and  1 1). 

Figures  3.94  through  3.127  contain  the  remaining  interfloat  travel  time  esti¬ 
mates  calculated  by  averaging  interpolated  reciprocal  path  measurements. 

3.1.2.  Sound  Speed  Estimates 

The  sound  speed  profile  for  the  upper  1000  meters  was  obtained  from  5 
expendable  bathythermograph  (XBT)  measurements  [13]  made  from  the  R/V  New 
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Horizon  near  the  deployment  site  at  various  times  during  the  Swallow  float  experiment 
and  historical  salinity  data  for  area  24,  3rd  quarter  of  the  year,  in  the  upper  1000  meters 
archived  by  the  National  Oceanographic  Data  Center  [13].  The  sound  speed  profile  for 
the  lower  3100  meters  was  obtained  from  a  conductivity-temperature-depth  CTD  station 
[14]  carried  out  on  July  11,  1989  near  the  center  of  the  float-triangle  where  all  the  freely 
drifting  floats  were  put  into  the  water.  Using  the  equation  given  by  Mackenzie  [15],  the 
composite  sound  speed  profile  was  calculated  and  plotted  in  Figure  3.128.  Harmonic 
mean  sound  speeds  for  each  interfloat  and  float  to  surface  path  were  also  calculated  and 
are  given  in  Appendix  1 . 

3.2.  Measurement  Statistics 

The  measurement  error  variance  '"timates  are  used  to  weight  the  measure¬ 
ments  so  that  greater  weight  is  given  to  measurement  with  a  smaller  error  variance.  Since 
travel  time  measurement  error  and  sound  speed  estimate  error  are  assumed  to  be  uncorre¬ 
lated,  the  measurement  error  variance  for  float  pair  ij  can  be  expressed  as  [4]: 

+  (4  y  (3.2) 

where  t,"  is  the  estimated  travel  time  between  float  i  and  j,  ol,  the  sound  speed  error  vari¬ 
ance  in  cij,  ci,  the  estimated  harmonic  sound  speed  between  float  i  and  j,  and  the  travel 
time  error  variance  between  float  i  and  j. 

3.2.1.  Travel  Time  Variances 

It  has  been  shown  [12]  that  the  variance  of  a  bottomed  float’s  surface  echo 
travel  time  measurements  can  be  used  as  an  estimate  of  the  variance  of  the  measurement 
error,  because  the  true  surface  echo  times  are  approximately  constant.  Subtracting  the 
means  from  the  surface  echo  travel  time  estimates  produces  estimates  of  the  travel  time 
estimate  errors.  The  variances  of  the  error  time  series  are  given  in  Table  3. 1 . 
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Table  3.1.  Estimated  Variance  of  Bottomed  Float 

Travel  Time,  July  1989  Experiment 

Float  Numbers 

Variance,  msec^ 

9 

9 

115 

10 

10 

no 

11 

11 

97 

9 

10 

78 

9 

11 

64 

10 

11 

142 

In  the  previous  section,  the  bottomed  float  surface  echo  travel  time  measure¬ 
ments  were  ruled  as  too  noisy  to  be  used  directly  in  estimating  the  float  positions,  and  the 
mode  of  the  travel  time  measurements  was  taken  as  the  travel  time  estimate.  The  vari¬ 
ance  of  the  error  in  the  constant  estimate  (mode)  is  expected  to  be  much  smaller  than  the 
variance  of  the  error  in  the  measurement.  Based  upon  the  predicted  maximum  vertical 
movement  of  a  tethered  bottom  float  and  wave  height  in  the  sea  surface  during  the 
deployment  (about  3  meters),  the  error  in  the  mode  is  taken  to  be  mean  zero  and  have 
variance  6  msec^.  By  the  same  token,  the  interfloat  travel  time  variance  is  taken  to  be  4 
msec2  [4]  based  upon  the  predicted  maximum  horizontal  movement  of  a  bottom  float. 
Table  3.2  summarizes  the  estimated  error  variances  of  the  bottomed  float  travel  time  esti¬ 
mates. 

The  surface  echo  travel  time  estimate  error  for  the  freely  drifting  floats  cannot 
be  estimated  directly.  Because  these  floats  were  deployed  shallower  than  the  bottom 
floats,  their  variances  are  expected  to  be  lower  than  the  values  in  Table  3.1  but  higher 
than  those  in  Table  3.2.  An  arbitrary,  yet  logical,  choice  of  error  variances  for  the  freely 
drifting  float  surface  echo  travel  time  estimates  are  listed  in  Table  3.3. 


12 


Table  3.2.  Estimated  Variance  of  Bottomed  Float 

Travel  Time  Based  on  Predicted  Float 

Movement,  July  1989  Experiment 

Float  Numbers 

Variance,  msec^ 

9 

9 

6 

10 

10 

6 

11 

11 

6 

9 

10 

4 

9 

11 

4 

10 

11 

4 

Table  3.3  Estimated  Variance  of  Freely  Drifting 

Float  Travel  Time,  July  1989  Experiment 

Float  Numbers 

Variance,  msec^ 

0 

0 

6 

1 

1 

10 

2 

2 

20 

3 

3 

30 

4 

4 

40 

5 

5 

50 

6 

6 

60 

7 

7 

70 

8 

8 

80 

Travel  time  error  variance  for  float  pairs  in  which  one  of  the  floats  is  not  sta¬ 
tionary  can  be  estimated  from  the  difference  between  reciprocal  path  travel  time  meas¬ 
urements  [12],  The  difference  in  clock  rates  causes  the  travel  time  difference  to  be  a 
linear  function  of  time.  The  difference  also  appears  to  contain  second  order  time  depen¬ 
dence  due  to  difference  in  clock  accelerations.  Subtracting  a  second  order  fit  from  the 
travel  time  difference  produces  an  estimate  of  the  error  in  the  travel  time  estimates.  Fig¬ 
ures  3.129  through  3.191  contain  interfloat  travel  time  differences  along  with  second 
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order  curves  fitted  to  the  differences  and  the  corresponding  error  estimate  time  series. 
The  variances  of  the  error  time  series  are  summarized  in  Appendix  2. 

3.2.2.  Sound  Speed  Variances 

It  has  been  shown  [2,4]  that  the  variance  of  the  sound  speed  estimate  at  a  par¬ 
ticular  depth  is  just  under  0.4  (meter/secp  in  the  North  Pacific  Ocean.  The  variance  of  the 
harmonic  mean  sound  speed  for  a  given  path  would  be  expected  to  be  less  than  that  in  the 
sound  speed  at  a  particular  depth  due  to  averaging  effect.  The  error  in  the  harmonic 
mean  sound  speed  is  estimated  to  be  0.1  (meter/secP  for  the  1989  experiment. 

3.3.  Initial  Estimates 

The  MPL  implementation  of  the  filters  requires  an  initial  estimate  of  the  float 
positions.  Additionally,  the  Kalman  filter  requires  an  initial  estimate  of  the  float  veloci¬ 
ties  and  the  initial  position  and  velocity  variance  estimates.  Before  estimating  the  initial 
positions,  we  must  first  define  a  coordinate  system.  The  coordinate  system  for  the  float 
localization  is  established  in  which  the  origin  lies  at  the  surface  directly  above  float  9. 
The  Z  axis  is  vertical,  positive  downwards  and  extends  through  float  9,  and  the  X  axis  is 
oriented  such  that  float  10  lies  in  the  XZ  plane.  The  positions  of  floats  9,  10  and  11  are 
taken  to  be  the  ship’s  position  when  the  floats  were  launched  into  the  water.  Figure 
3.192  shows  a  plane  view  of  the  coordinate  system.  The  locations  of  the  freely  drifting 
floats  will  be  estimated  relative  to  these  axes. 

3.3.1.  Initial  Position  Estimates 

The  time  for  the  initial  float  position  chosen  for  float  localization  is  record 
1003,  00.00  July  9  1989,  at  which  7  of  the  9  freely  drifting  floats  had  reached  equilibrium 
depth.  The  MPL  implementation  of  the  least  squares  filter  requires  a  rough  initial  posi¬ 
tion  estimate  to  bootstrap  the  filter.  The  rough  initial  positions  are  taken  to  be  the  ship’s 
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position  relative  to  the  coordinate  system  when  the  floats  were  launched  into  the  water 
using  depths  taken  from  Figures  3.1  through  3.12.  Table  3.6  lists  the  rough  initial  posi¬ 
tion  estimate. 


Table  3.6.  Rough  Initial  Position  Estimate 

Record  I003.00:00July  9  1989 

Roat 

X 

Y 

Z 

9 

0 

0 

4050 

10 

6150 

0 

4050 

11 

3075 

5500 

4050 

0 

3075 

2750 

570 

1 

3075 

2750 

990 

2 

3075 

2750 

1380 

3 

3075 

2750 

1670 

4 

3075 

2750 

2080 

5 

3075 

2750 

2490 

6 

3075 

2750 

2840 

7 

3075 

2750 

3350 

8 

3075 

2750 

3650 

Using  the  rough  initial  position  estimates  in  Table  3.6,  the  travel  time  estimates 
for  record  1003,  the  travel  time  variance  estimates  given  in  Tables  3.2  and  3.3  and 
Appendix  2,  and  the  sound  speeds  and  sound  speed  variances  given  in  Appendix  1,  the 
least  squares  filter  produced  the  position  estimate  given  in  Table  3.7  It  has  a  root  mean 
squared  (rms)  residual  of  2.67  meters.  Since  an  initial  position  estimate  is  determined  to 
be  satisfactory  when  the  least  squares  filter  converges  to  a  position  estimate  which  results 
in  an  rms  residual  of  less  than  7.5  meters  [4],  the  position  estimate  at  record  1003  is 
taken  as  a  good  initial  position  estimate  and  will  be  used  by  the  Kalman  filter. 
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Table  3.7.  Initial  Position  Estimate  Produced 
by  the  Least  Squares  Filter  for  Record  1003, 
00:00,  July  9  1989 

Float 

X 

Y 

Z 

9 

0 

0 

4039 

10 

6131 

0 

4051 

11 

3103 

5518 

4041 

0 

3779 

-1456 

574 

1 

4887 

82 

999 

2 

3885 

162 

1380 

3 

4425 

969 

1674 

4 

4110 

1398 

2081 

5 

4092 

1482 

2489 

6 

3870 

1754 

2841 

7 

2875 

2085 

3353 

8 

2824 

1871 

3650 

3.3.2.  Initial  Velocity  Estimates 

Using  the  good  initial  position  estimates,  the  least  squares  filter  was  again  run 
for  the  first  13  records  (9  minute  period)  from  records  1003  to  1015.  The  position  esti¬ 
mates  produced  by  the  filter  were  used  to  calculate  the  position  change  rates 
meters/second.  Twelve  such  successive  position  change  rates  were  then  averaged  to  pro¬ 
duce  the  initial  velocity  estimates.  Table  3.8  lists  the  velocity  estimates  for  record  1003, 
00:00,  9  July  1989. 

3.3.3.  Initial  Position  and  Velocity  Variance  Estimates 

The  estimate  of  the  variances  in  the  initial  positions  and  velocities  is  also 
needed  by  the  Kalman  filter.  The  estimate  suggested  by  [4]  will  be  used  and  is  listed  in 
Table  3.9. 
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Table  3.8.  Initial  Velocity  Estimate  (in  meters/second). 

Record  1(X)3,  (X):(X),  July  1989  Experiment 

Float 

X 

Y 

z 

9 

0 

0 

0 

10 

0 

0 

0 

11 

0 

0 

0 

0 

0.061 

-0.040 

-0.004 

1 

0.025 

-0.058 

0.005 

2 

-0.015 

-0.040 

0.001 

3 

0.028 

-0.027 

0.006 

4 

0.003 

0.001 

0.013 

5 

-0.003 

0.001 

0.015 

6 

0.002 

0.017 

0.027 

7 

-0.010 

0.021 

0.005 

8 

-0.036 

-0.003 

0.025 

Table  3.9.  Estimate  of  Initial  Position  and  Velocity  Variance 
Used  by  the  Kalman  Filter,  July  1989  Experiment 

Position,  meters^ 

Velocity,  (m/sec)2 

Float 

X 

Y 

z 

X 

Y 

z 

9 

0 

0 

0.01 

0 

0 

0 

10 

0.01 

0 

0.01 

0 

0 

0 

11 

0.01 

0.01 

0.01 

0 

0 

0 

0 

0.01 

0.01 

0.01 

0.0001 

0.0001 

0.00001 

1 

0.01 

0.01 

0.01 

0.0001 

0.0001 

0.00001 

2 

0.01 

0.01 

0.01 

0.0001 

0.0001 

0.00001 

3 

0.0 1 

0.0 1 

0.01 

0.0001 

0.0001 

0.00001 

4 

0.01 

0.01 

0.01 

0.0001 

0.0001 

0.00001 

5 

0.01 

0.01 

0.01 

0.0001 

0.0001 

0.00001 

6 

0.01 

0.01 

0.01 

0.0001 

0.0001 

0.00001 

7 

0.01 

0.01 

0.01 

0.0001 

0.0001 

0.00001 

8 

0.01 

0.01 

0.01 

0.0001 

0.0001 

0.00001 

3.4.  Process  Noise 

The  last  piece  of  information  required  by  the  Kalman  filter  is  the  float  accelera¬ 
tion  variance,  also  known  as  the  process  noise,  described  in  equation  (2.6).  The  process 
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noise  is  a  parameter  used  by  the  Kalman  filter  to  determine  how  much  weight  to  give  to 
its  own  track  of  the  float  relative  to  the  measurements.  When  the  model  (i.e.  filter)  pro¬ 
cess  noise  matches  the  true  process  noise  experienced  by  the  float  during  the  experiment, 
the  Kalman  filter  is  performing  optimally  and  will  have  more  estimates  falling  inside  the 
predefined  confidence  interval  [16].  However,  there  is  no  easy  way  of  estimating  the  true 
process  noise.  One  approach  suggested  by  [4]  is  to  run  the  filter  several  times,  each  time 
with  a  different  process  noise  value  and  the  one  which  minimizes  the  innovation  power 
will  be  selected  as  the  candidate.  We  will  use  this  matched  processing  technique  with 
values  ranging  from  l0->2  to  5  x  lO-s  (meter/second2)2  in  search  of  the  true  process  noise. 
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4.  Filter  tuning 


4.1.  Least  Squares  Filter 

Using  the  good  initial  positions  given  in  Table  3.7,  the  sound  speeds  and  sound 
speed  variances  given  in  Appendix  1,  the  travel  time  estimates  for  records  1003  to  2120, 
and  the  travel  time  variances  given  in  Tables  3.2  and  3.3  and  Appendix  2,  the  least 
squares  filter  produced  position  estimates  with  an  average  rms  of  2.7  meters.  The  residual 
is  low  and  indicates  that  the  position  estimate  is  acceptable.  However,  a  close  inspection 
of  the  residual  sequences  Z,  -  /i(X,)]  reveals  that: 

I 


E  Z  (iirtct  path  b£tw€e/i  9  A.  10  ~  ^  path  bativaan  9  4  /O  )  ~  1*3  ITlCtCrS 

E  Z  diract  path  bttwaen  9  A  U  “  ^^Xdir*ctpathb€twatii9  A  //  )  “  3.0  ITlCtCrS 
E  Zdxrtct  pathbatwtan  10  A  //  ^  diract  path  /O  dr  /  /  )  “  1*2  rnCtCFS 

£  ^ flo<u  9  surfaca  acho  “  hijC float  9  aurfaca  echo  )  —1.5  mCtCfS 

E  Z  float  10  surfaca  acho  ^  float  W  suffaca  acho')  ~  4,9  ITlCtCrS 
E  Z float  II  surfaca  acho  ~  h  (x float  It  surface  acho  )  “  0.6  ITIClCrS 

By  definition,  residuals  in  the  measurement  model  are  to  be  mean  zero: 


E  V,  =e[z, =0  (4.7) 

The  bias  contained  in  the  residual  sequences  (4.1)  -  (4.6)  is  an  indication  that  the  mode 
values  selected  as  the  travel  time  estimates  may  not  be  the  optimal  choice  for  the  July 
1989  data  set.  In  order  to  improve  the  filter  performance,  trends  in  the  residual 
sequences  must  be  removed.  An  iterative  procedure  was  developed  to  achieve  this  and  is 


described  as  follows: 
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1. 

2. 

3. 

4. 


Compute  At,/  = 


2ii  -  hixPj) 


where  i  =  9,10,1 1  and  j  =  9,10,1 1. 


Update  T,/  =  t,/  -  At./  . 

Rerun  the  least  squares  filter. 


If  all 


E  2i/  -  h(x!j) 


<  1  meter  then  stop;  otherwise  go  to  step  1 . 


The  procedure  was  applied  to  the  July  1989  data  set.  After  10  iterations,  all  6  expected 
residuals  converged  to  within  1  meter,  and  the  least  squares  filter  produced  the  position 
estimate  with  an  average  rms  of  2.3  meters.  The  residual  is  lower  than  2.7  meters  and 
indicates  that  the  position  estimate  for  records  1003  to  2120  is  closer  to  the  true  float 
positions.  The  constant  value  estimates  used  in  the  10th  run  listed  in  Table  4.1  will 
replace  the  mode  values  as  the  new  travel  time  estimates  and  will  be  used  by  the  Kalman 
filter. 


Table  4.1.  Bottomed  Float  Surface  Echo  and 

Direct  Path  Between  Bottomed  Floats  Travel 

Time  Estimate,  July  1989  Experiment 

Float  Path 

Travel  Time  Estimate  (seconds) 

9-10 

4.028  (was  4.025) 

9-11 

4.152  (was  4.157) 

10-11 

4.130  (was  4.132) 

9-9 

5.401  (was  5.402) 

10-10 

5.409  (was  5.422) 

11-11 

5.405  (was  5.404) 

4.2.  Kalman  Filter 

The  measurements,  sound  speed  and  variance  estimates,  and  the  initial  esti¬ 
mates  were  applied  to  the  Kalman  filter.  Using  process  noise  values  ranging  from  10-';  to 
5xio-s  (meters/second2)2,  the  Kalman  filter  was  run  24  times.  Table  4.2  lists  the  mean  inno- 
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vation  power  resulting  from  each  run. 


Table  4.2.  Mean  Innovation  Power  Produced 

by  the  Kalman  Filter.  July  1989  Experiment 

Process  Noise 

Mean  Innovation  Magnitude 

(meters/second2)2 

(meters) 

1x10-12 

7.0671 

2x10-12 

6.0969 

3x10-12 

5.6302 

5x10-12 

5.1290 

7x10-12 

4.8450 

1x10-11 

4.5793 

2x10-11 

4.1538 

3x10-11 

3.9536 

5x10-11 

3.7460 

7x10-11 

3.6335 

1x10-10 

3.5331 

2x10-10 

3.3866 

3x10-10 

3.3269 

5x10-10 

3.2753 

7x10-10 

3.2546 

1x10-9 

3.2435 

2x10-9 

3.2539 

3x10-9 

3.2808 

5x10^ 

3.3401 

7x10-9 

3.3974 

lxl0-« 

3.4770 

2xl0-« 

3.6933 

3xl0-» 

3.9648 

5x10-* 

6.8216 

It  has  been  shown  using  simulations  [4]  that  the  process  noise  which  minimizes  the 
power  in  the  innovations  sequence  produces  a  position  estimate  with  the  smallest  true 
rms  error.  The  minimum  is  seen  to  occur  at  lO-’  which  is  the  same  found  for  the  1987 
experiment  [4].  Thus,  process  noise  for  the  July  1989  experiment  is  estimated  to  be 


lO-’  (m/s02. 
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5.  Discussion  of  Localization  Results 

In  this  section,  the  outputs  from  the  localization  filters  are  discussed.  Of 
interest  are: 

•  RMS  position  error  estimate 

•  RMS  residual/innovation 

•  Float  position  estimate 

•  Distance  difference  between  filters 

•  Float  depth  variation 

•  Float  speed  estimate 

The  curves  shown  in  Figure  5.1  are  the  mean  innovation  power  1 1 1,  the  Kalman 
filter’s  estimate  of  the  rms  position  error,  and  the  least  squares  filter’s  estimate  of  the  rms 
position  error.  The  vertical  axis  is  rms  error  or  mean  innovation  magnitude  in  meters. 
The  horizontal  axis  is  the  estimate  of  process  noise  given  to  the  Kalman  filter  as  a  param¬ 
eter.  Simulations  [2]  have  shown  that  the  true  rms  error  in  the  Kalman  filter’s  position 
estimate  was  always  bracketed  by  the  mean  magnitude  of  the  innovation  and  the  filter’s 
estimate  of  the  rms  error.  Therefore,  the  rms  error  in  the  Kalman  filter  estimate  is 
thought  to  be  less  than  3.3  meters.  The  least  squares  filter’s  estimate  of  rms  float  position 
error  is  2.3  meters.  Simulation  results  have  indicated  that  the  true  rms  error  in  the  least 
squares  filter’s  position  estimate  is  no  more  than  about  twice  the  estimated  error. 
Accordingly,  the  rms  error  in  the  least  squares  filter  estimate  is  thought  to  be  less  than  4.6 
meters. 

The  rms  residual  for  the  least  squares  estimates  and  the  rms  innovation  for  the 
Kalman  filter  estimates  are  shown  as  a  function  of  measurement  (Swallow  float  record) 
number  in  Figure  5.2.  Both  rms  residual  and  innovation  contain  large  spikes  around 
records  1600  and  1700.  These  spikes  are  due  to  missing  surface  echo  travel  time  meas- 
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urements  in  floats  6  and  8’s  data  and  the  values  determined  by  the  adaptive  linear  predic¬ 
tor  must  still  contain  large  errors. 

Figures  5.3a  and  5.3b  show  the  filters’  float  horizontal  position  estimates 
between  records  1003  and  2120.  Both  position  estimates  indicate  that  the  freely  drifting 
floats  dispersed  away  from  the  center  of  the  float  triangle  with  floats  0  and  1  moving  to 
the  northwest,  floats  2  and  3  to  the  west,  float  4  to  the  southwest,  and  floats  5,  6,  7  and  8 
to  the  southeast.  The  drifting  pattern  was  probably  due  to  the  complex  water  movement 
(including  eddies,  wind-driven  surface  currents  and  the  California  Current)  near  the 
experiment  site.  The  filters’  float  depth  estimates  are  also  plotted  in  Fig  5.4a  and  5.4b. 

The  distance  between  the  two  filters’  estimates  is  plotted  for  each  float  as  a 
function  of  measurement  number  in  Figure  5.5.  The  curves  are,  in  descending  order,  from 
floats  0  through  1 1.  The  distance  between  the  two  position  estimates  is  relatively  small 
for  all  floats  with  an  rms  difference  less  than  3  meters  for  most  of  the  experiment. 

Figures  5.6  through  5.8  show  the  enlarged  version  of  the  Kalman  filter’s  float 
depth  estimates  as  a  function  of  measurement  number.  The  floats  oscillate  irregularly 
about  slowly  varying  means  with  periods  on  the  order  of  30  to  90  minutes  which 
corresponds  roughly  to  internal  wave  periods.  The  slowly  increasing  depth  trend,  except 
for  float  0,  is  probably  caused  by  gradual  compression  of  the  float  as  it  gets  colder. 

The  float  speed  estimates  derived  from  the  Kalman  filter’s  X,  Y,  and  Z  velocity 
component  estimates  are  plotted  in  Figures  5.9  through  5.1 1.  The  float  speeds  appear  to 
be  rather  unsteady  (i.e.,  high  process  noise)  throughout  the  experiment.  The  large  scale 
speed  oscillation  with  periods  on  the  order  of  4  to  12  hours,  depending  on  float  equili¬ 
brium  depth,  is  thought  to  be  caused  by  the  tidal  currents.  Average  float  speeds  for  the 
freely  drifting  floats  are  listed  in  Table  5.1. 


6,  Summary 

In  this  report,  two  localization  methods,  generalized  least  squares  filter  and 
Kalman  filter,  were  reviewed  from  the  performance  criteria  point  of  view.  The  inputs  to 
both  filters  using  data  from  the  July  1989  experiment  were  described  in  detail.  The  pro¬ 
cedures  for  tuning  the  filters  were  described  and  the  results  from  both  filters  were  dis¬ 
cussed  and  compared.  The  rms  position  error  is  estimated  to  be  less  than  3.3  m  for  the 
Kalman  filter  and  less  than  4.6  m  for  the  least  squares  filter.  Due  to  the  high  process 
noise  experienced  by  the  floats  during  the  experiment  and  the  quasi  vertical  line  array 
deployment  geometry  [3],  the  two  filters  perform  comparably.  The  two  position  esti¬ 
mates  have  an  rms  difference  of  less  than  3  meters  for  most  of  the  experiment.  Both 
methods  appear  to  be  capable  of  estimating  float  positions  to  within  the  desired  accuracy 
of  one-tenth  of  a  wavelength  at  the  highest  frequency  of  interest  20  Hz  (7.5  m)  in  order  to 
effectively  beamform  the  VLF  acoustic  data. 
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Appendix  1  -  Harmonic  Mean  Sound  Speed  Estimates 


Harmonic  Mean  Sound  Speed  Estimates 
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Harmonic  Mean  Sound  Speed  Estimates 
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Sound  Speed,  tn/msec 

Error  Variance,  (m/msec)2 
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Figure  3.2 
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Figure  3.3 


Pulse  Leading  Edges.  Fioa-t  3  listening  to  itself,  July  1969 
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Figure  3.4 
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Figure  3.5 


Pulse  Leading  Edges.  Floai  6  listening  to  itself,  July  1989 
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Figure  3.7 


Pulse  Leading  Edges.  Flcni  7  listening  to  itself,  July  1989 
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Pulse  Leading  Edges.  Floa't  6  listening  to  itself,  July  1989 
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Figure  3.9 


Pulse  Leading  Edges.  Floa't  9  listening  to  itself,  July  1989 


Figure  3.10 


Pulse  Leading  Edges.  Floa-t  10  listening  to  itself,  July  1989 
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Figure  3.11 


Pulse  Lecxding  Edges.  Flca-t  11  listening  to  itself,  July  1969 


Figure  3.12 
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Pulse  Leading  Edges.  Flcai  9  listening  to  Float  10,  July  1989 
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Pulse  Leading  Edges.  Float  10  listening  to  Float  9,  July  1989 
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Figure  3.13 
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Pulse  Leading  Edges.  Floai  9  listening  to  Float  11,  July  1989 
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Pulse  Leading  Edges.  Float  11  listening  to  Float  9,  July  1989 
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Figure  3.14 


Pulse  Leading  Edges.  Floa't  9  listening  to  Float  L,  July  1989 
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Pulse  Leading  Edges.  Float  1  listening  to  Float  9,  July  1989 
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Figure  3.16 


Pulse  Leading  Edges.  Float  9  listening  to  Float  2,  July  1969 
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Pulse  Leading  Edges.  Float  2  listening  to  Floot  9,  July  1989 

I - 1 - 1 - 1 - 1 - 1 - 1 - 1 - y 


I - 1 - 1 - 1 - \ - 1 - 1 - 1 - 1— J 

600  800  1000  1200  1400  1600  1800  2000  2200 

record  nmiber 


Figure  3.17 


Pulse  Leading  Edges.  Float  9  listening  to  Floot  3,  July  1989 


Pulse  Leading  Edges.  Float  3  listening  to  Float  9,  July  1989 


Figure  3.18 
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Pulse  Leading  Edges.  Flcxit  9  listening  to  Float  4,  July  1989 
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Figure  3.19 
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Pulse  Leading  Edges.  Floal  9  listening  to  Float  6.  July  1989 
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Pulse  Leading  Edges.  Float  6  listening  to  Float  9,  July  1989 
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Figure  3.21 


Pulse  Leading  Edges.  Floai  ^  listening  to  Float  7,  July  1989 


Pulse  Lending  Edges.  Float  7  listening  to  Float  9,  July  1989 
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Figure  3.22 


Pulse  Leading  Edges.  Flcxi't  10  listening  to  Float  11,  July  1989 
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Pulse  Leading  Edges.  Float  11  listening  to  Floot  10,  July  1989 
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Figure  3.24 
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Pulse  Leading  Edges.  Floai  10  listening  to  Float  0,  July  1989 
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Pulse  Leading  Edges.  Float  0  listening  to  Float  10,  July  1989 


Figure  3.25 
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Pulse  Leading  Edges.  Flcni  10  listening  to  Float  1,  July  1989 
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Pulse  Leading  Edges.  Float  1  listening  to  Float  10,  July  1959 
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Figure  3.26 


Pulse  Leading  Edges.  Float  10  listening  to  Float  3,  July  1989 


Pulse  Leading  Edges.  Float  3  listening  to  Float  10.  July  1989 


Figure  3.28 
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Pulse  Leading  Edges.  Float  10  listening  to  Float  4,  July  1989 


Pulse  Leoding  Edges.  Float  4  listening  to  Float  10,  July  1989 
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Figure  3.29 
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Pulse  Leading  Edges.  Float  10  listening  to  Float  5,  July  1989 
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Figure  3.30 
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Pulse  Leading  Edges.  Float  10  listening  to  Float  6,  July  1989 


Figure  3.31 
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Pulse  Leading  Edges.  Float  10  listen-ing  to  Float  7,  July  1989 
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Pulse  Leading  Edges.  Float  7  listening  to  Floot  10,  July  1989 
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Figure  3.32 


Pulse  Leading  Edges.  Float  10  listeaing  to  Float  8.  July  1989 
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Pulse  Leading  Edges.  Float  8  listening  to  Float  10,  July  1989 


Figure  3.33 
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Pulse  Leading  Edges.  Floai  11  listening  to  Float  0,  July  1969 


Pulse  Leading  Edges.  Float  0  listening  to  Float  11,  July  1989 
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Figure  3.34 


Pulse  Leading  Edges.  Float  11  listening  to  Float  1,  July  1989 


Pulse  Leading  Edges.  Float  1  listening  to  Floot  11,  July  1989 
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Figure  3.35 
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Pulse  Leading  Edges.  Flcat  11  listening  to  Float  2,  July  1989 
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Pulse  Leading  Edges.  Float  2  listening  to  Float  11,  July  1989 
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Figure  3.36 
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Pulse  Leading  Edges.  Floa*  11  listening  to  Float  3,  July  1989 


Pulse  Leading  Edges.  Float  3  listening  to  Float  11.  July  1989 
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Figure  3.37 
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Pulse  Leading  Edges.  Float  11  listening  to  Float  4,  July  1989 
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Pulse  Leading  Edges.  Float  4  listening  to  Float  11,  July  1989 
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Pulse  Leading  Edges.  Float  11  listening  to  Float  5,  July  1989 


600  800  1000  1280  1400  1600  1800  2000  2208 


record  number 


Pulse  Leading  Edges.  Flcat  5  listening  to  Float  11,  July  1989 


record  nu.’^cer 


Figure  3.39 


Pulse  Leading  Edges.  Float  11  listening  to  Float  6,  July  1989 


Pulse  Leoding  Edges.  Float  6  listening  to  Float  11.  July  1989 


Figure  3.40 
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Pulse  Leading  Edges.  Float  11  listenir-  to  t- loot  7,  July  1959 
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Pulse  Leading  Edges.  Fioct  11  listening  to  Float  8,  July  1989 


Pulse  Leoding  Edges.  Float  8  listening  to  Float  11.  July  1989 
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Figure  3.42 
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Pulse  Lecxjfng  Edges.  Flca't  0  listening  to  Float  1,  July  1989 
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Pulse  Leading  Edges.  Float  1  listening  to  Float  0,  July  1989 


Figure  3.43 
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Pulse  Leading  Edges.  Floa-t  0  listening  to  Float  2,  July  1969 
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Figure  3.44 
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Pulse  Leading  Edges.  Float  0  listening  to  Float  3,  July  1989 
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Pulse  Leading  Edges.  Float  3  listening  to  Float  0,  July  1989 
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Figure  3.45 
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Figure  3.46 


Pulse  Leading  Edges.  Floa-t  0  listening  to  Float  5,  July  1989 
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Pulse  Leading  Edges.  Floai  0  listening  to  Float  6,  July  1969 
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Pulse  Leading  Edges.  Float  6  listening  to  Float  0,  July  1969 
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Figure  3.48 
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Pulse  Leading  Edges.  Float  0  listening  to  Float  7,  July  1969 
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Pulse  Leading  Edges.  Float  7  listening  to  Float  0,  July  1989 
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Figure  3.49 
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Pulse  Leodfng  Edges.  Float  1  listening  to  Float  2.  July  1989 
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Pulse  Leading  Edges.  Float  2  listening  to  Float  1,  July  1989 
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Figure  3.51 
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Pulse  Leading  Edges.  Floai  1  listening  to  Float  3,  July  1989 
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Pulse  Leading  Edges.  Float  3  listening  to  Float  1,  July  1989 
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Figure  3.52 


Pulse  Leading  Edges.  Floa't  1  listening  to  Float  4,  Jul\/  1989 
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Figure  3.53 


Pulse  Leading  Edges.  Floai  1  listening  to  Float  5,  July  1989 
I - 1 - 1 - 1 - 1 - 1 - 1 - 1 - h 


I - 1 - 1 - 1 - 1 - 1 - 1 - 1 - h-J 

600  800  1000  1200  1400  1600  1800  2000  2200 

record  runber 

Pulse  Leading  Edges.  Float  5  listening  to  Floot  1,  July  1989 
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Figure  3.54 
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Pulse  Leading  Edges.  Floot  1  listening  to  Float  6.  July  1989 
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Pulse  Leading  Edges.  Float  6  listening  to  Float  1,  July  1989 
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Pulse  Leading  Edges.  Float  1  listening  to  Float  7,  July  1989 
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Pulse  Leading  Edges.  Float  7  listening  to  Float  1,  July  1989 
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Figure  3.56 


Pulse  Leading  Edges.  Flood  1  listening  to  Float  8,  July  1989 
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Pulse  Leading  Edges.  Float  8  listening  to  Float  1,  July  1989 
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Figure  3.57 
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Figure  3.58 


Pulse  Leading  Edges.  Floa-t  2  listening  to  Float  4.  July  1989 


H - 1 - 1 - 1 - 1 - 1 - 1 - h 


I - 1 - 1 - 1 - 1 - 1 - ^ - 1 - h-* 

600  800  1000  1200  1400  1600  1800  2000  2200 

record  minber 


Pulse  Leading  Edges.  Float  4  listening  to  Floot  2,  July  1989 


Figure  3.59 


Pulse  Leading  Edges.  Floa-t  2  lis-tening  to  Float  5,  July  1989 
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Pulse  Leading  Edges.  Float  5  listening  to  Float  2,  July  1989 


record  runber 


Figure  3.60 
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Pulse  Lecxiing  Edges.  Float  2  listening  to  Float  6,  July  1969 
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Pulse  Leading  Edges.  Float  6  listening  to  Float  2,  July  1989 
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Figure  3.61 
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Pulse  Leading  Edges.  Flco't  2  listening  to  Float  7,  July  1989 
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Pulse  Leading  Edges,  Float  7  listening  to  Float  2,  July  1989 
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Figure  3.62 
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Pulse  LecxJing  Edges.  Flcn-t  2  listening  to  Float  9,  July  1999 
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Pulse  Leading  Edges.  Float  9  listening  to  Float  2,  July  1999 
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Figure  3.63 
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Pulse  Leading  Edges.  Flcat  3  lis-tening  to  Float  4,  July  1989 


Pulse  Leading  Edges.  Float  4  listening  to  Floot  3,  July  1989 


Figure  3.64 


Pulse  Leading  Edges.  Float  3  listening  to  Float  5,  July  1989 
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Pulse  Leading  Edges.  Float  5  listening  to  Float  3,  July  1989 
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Figure  3.65 
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Pulse  Leading  Edges.  Floa-t  3  listening  to  Float  6,  July  1989 


Pulse  Leading  Edges.  Float  6  listening  to  Float  3,  July  1989 
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Pulse  Leading  Edges.  Float  7  listening  to  Float  3,  July  1989 
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Figure  3.67 
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Pulse  Leading  Edges.  Floa-t  3  listening  to  Float  8,  July  1989 
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Pulse  Leading  Edges.  Float  8  listening  to  Float  3,  July  1969 
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Figure  3.68 
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Pulse  Leading  Edges.  Floai  4  listening  to  Float  5,  July  1989 


600  800  1000  1200  1400  1600  1800  2000  2200 


record  ntrober 

Pulse  Leading  Edges.  Float  5  listening  to  Float  4,  July  1989 
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Figure  3.69 
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Pulse  Lecxding  Edges,  FloQ't  4  listening  to  Float  6,  July  1969 
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Pulse  Leading  Edges.  Float  6  listening  to  Floot  4.  July  1989 
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Pulse  Leading  Edges.  Flcat  4  listening  to  Float  7,  July  1989 
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Pulse  Leading  Edges.  Float  7  listening  to  Float  4,  July  1989 

S 

I 

|i 

I 

i 

600  800  1000  1200  1400  1600  1800  2m3  2200 

record  number 


Figure  3.71 
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Pulse  Leading  Edges.  Floai  4  lis-tening  to  Float  6,  July  1969 
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Pulse  Leading  Edges.  Float  8  listening  to  Float  4,  July  1989 
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Figure  3.72 
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Pulse  Leading  Edges.  Float  5  listening  to  Floot  6,  Juli/  1959 
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Pulse  Leading  Edges.  Float  6  listening  to  Float  5,  July  1989 
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Figure  3.73 


Pulse  Lecxiing  Edges.  Flcai  5  listening  to  Float  7,  July  1989 
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Figure  3.74 
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Pulse  Leading  Edges.  Flooi  5  listening  to  Float  6,  July  1989 
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Pulse  Leading  Edges.  Float  8  listening  to  Float  5,  July  1989 
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Figure  3.75 
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Pulse  Leading  Edges.  Floai  6  listening  to  Float  7,  July  1989 
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Pulse  Leading  Edges.  Float  7  listening  to  Float  6,  July  1989 
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Figure  3.76 


Pulse  Leading  Edges.  Float  6  listening  to  Float  8,  July  1989 
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Pulse  Leading  Edges.  Float  8  listening  to  Float  6,  July  1989 
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Figure  3.77 
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Pulse  Leading  Edges.  Floai  7  listening  to  Float  8,  July  1989 
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Pulse  Leading  Edges.  Float  8  listening  to  Float  7,  July  1989 


Figure  3.78 
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Surface  echo  -travel  time  estimate,  float  0,  July  1989 
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Figure  3.79 
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Surface  echo  -travel  time  estimate,  float  2,  July  1989 
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Surface  echo  travel  time  estimate,  float  3,  July  1989 
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Surface  echo  travel  'time  estimate,  float  6,  July  1989 
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Figure  3.85 
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Surface  echo  travel  time  estimate,  float  7,  July  1989 
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Figure  3.86 
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Surface  echo  travel  time  estimate,  float  11,  July  1989 
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Float  11  histogram  of  surface  echo  travel  time  minus  the  mean  5411.8, 
July  1989  exp^iment 
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Histogram  of  travel  time  minus  the  mean  4017.7,  floats  9  and  10 
July  1989  experiment 
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Histogram  of  travel  time  minus  the  mean  4030.2.  floats  9  ond  10 
July  1989  experiment 


Figure  3.91c 
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Travel  time  estimate,  floats  9  and  11,  July  1989 

Direct  path,  calculated  using  smoothed  surface  echo  travel  times 
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Histogram  of  travel  time  minus  the  mean  4159.4,  floats  9  and  11 
July  1989  experiment 
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Trayel  time  estimate,  floats  10  and  11,  July  1989 
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Figure  3.93a 


Travel  time  estimate,  floats  10  and  11.  July  1989 

Diroct  path,  calculated  using  actual  sirtacs  gox3  travel  times 
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Histogrcm  of  travel  time  minus  the  mean  4L29.6.  floats  10  end  11 
July  1989  experiment 
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Histogram  of  travel  time  minus  the  mean  4137.3,  floats  10  and  11 
July  1989  experiment 
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Travel  time  estimate,  floats  9  and  0,  July  1989 
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Travel  time  e<"'‘  ■  ,,  flaats  9  and  1.  July  1989 
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Travel  time  estimate,  floats  9  and  2,  July  1989 
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Travel  time  astimate,  floats  9  and  5.  July  1989 
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Travel  time  estimate,  flaats  9  and  6,  July  1989 


msec 


2543. 

2362. 

2180. 

1998. 

1817. 

1635. 


Travel  time  esiinrate,  floats  10  ond  0,  July  1989 


Travel  time  estimate,  flaats  10  and  1,  July  1989 

I - 1 - 1 - h 


2520. 


2450. 

2380. 

oi 

1/1 

E  2310. 

2240. 

2170. 


I - 1 - 1 - \ - h  -  - 

600  800  1000  1200  1400  1600  1 

record  HLinber 


Figure  3.99 
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Travel  time  estimate,  flaats  11  and  4,  July  1989 
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Figure  3. 109 
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Travel  time  estimate,  floats  4  and  8,  July  1969 
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Figure  3. 123 
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Travel  time  estimate,  floats  5  and  7,  July  1989 
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Travel  time  estimate,  floats  5  and  8.  July  1989 
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Figure  3. 125 
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Travel  time  esiitnate,  floats  7  and  8,  July  1989 
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Travel  time  difference  minus  second  order  fit,  floats  9  and  0 
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Figure  3. 129 
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Travel  time  difference  minus  second  order  fit,  floats  9  and  5 
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Trovel  time  difference  minus  second  order  fit,  floats  9  and  7 
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Travel  time  difference  minus  second  order  fit.  floats  10  and  7 
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Travel  time  difference  minus  second  order  fit,  floats  11  and  0 
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Travel  time  difference  mitius  second  order  fit,  floats  11  and  5 
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Travel  time  difference  minus  second  order  fit,  floots  11  and  7 
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Figure  3. 155 
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Travel  time  difference  minus  second  order  fit,  floats  0  and  1 
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Figure  3. 156 
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Travel  time  difference  minus  second  order  fit,  floats  0  and  3 
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Figure  3. 158 
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Travel  time  difference  minus  second  order  fit.  floots  0  and  5 
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Figure  3. 160 
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Travel  time  difference  minus  second  order  fit,  floats  1  and  4 
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Figure  3. 166 
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Travel  time  difference  minus  second  order  fit,  floots  1  and  7 
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Travel  time  difference  minus  second  order  fit,  floats  2  and  5 
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Travel  time  difference  minus  second  order  fit,  floats  2  ond  6 
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Fieiire  3. 174 
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Travel  time  difference  minus  second  order  fit,  floats  2  ond  9 
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Figure  3. 176 


Travel  time  difference,  floats  3  and  8. Juh/  1999.  md  a 

second  order  fit)  coefficients  oroi  -0.02307  -23.57817  and  35.45761 
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Travel  time  difference  minus  second  order  fit,  floats  3  and  8 
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Travel  time  difference  minus  second  order  fit,  floats  4  and  5 
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Travel  time  difference,  floats  4  and  6.  July  19Q2»  end  a 
second  order  fit|  coefficients  orsi  (none)  -2o.25o%  end  -64.93631 
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Travel  time  difference  minus  second  order  fit,  floats  4  and  6 
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Figure  3. 183 
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Travel  time  difference,  floats  5  and  l‘^9,_md  a 

sectxid  order  fit)  coefficients  arsi  -0.0576o  -35. 127W  and  20.22412 
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Travel  time  difference  minus  second  order  fit,  fleets  5  ond  7 
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Fiiiure  3. 187 
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Figure  3. 189 
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Plane  view  of  coordinate  system  used  for  float  localization 


Figure  3. 192 


Mean  Irrovation  or  RMS  Error  (meter) 


Results  of  the  Kalman  and  Least  Squares  Filters,  July  1989  Experiment 
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Figure  5.4b 
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Floo^  d  dapth  asitmates,  Kalman  FlHar,  July  1989 
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Figure  5.6 


rioai  6  depth  eetimatee,  Kalman  Filter,  July  1939 


Float  7  depth  eetimatee,  Kalman  Flltw^,  July  1939 


Float  8  depth  eetimatee,  Kalman  Filter,  July  1939 


